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Waste Savings Analysis 
Exemplary Newspaper Application 

1900 


Customer Information 


Description 


Value 


1901 Circulation 


50000 


1902 Width (inches) 


14.75 


1903 Length (inches) 


11.376 


1 904 Paper Pages Across Web Width 


2 


1 906 Paoer Paaes Around Circumference 


4 


1 906 # of Paper Pages/Cylinder Revolution 


8 


1 907 Newsprint Cost (US$/Pound) 


US$ 0.25 


1 908 Paper Weight (Pounds/Ream) 


30 


1909 Ink Cost (Percent of Paper Cost) 


10% 


1910 Mal<e Ready Count / Day On Edition 


3 


1911 Days per Week 


6 


1912 Specials per Week 


5 


1913 Number of Paper Pages per Newspaper 


36 


1914 Savings (Papers / Make Ready) 


1000 


Calculations 


1915 Number of Webs 


4.5 


1916 Sauare Feet of Paper per Newspaper 


41.95 


1917 Weight/Paper (Pounds) 


0.39 


1918 Paper Cost per Newspaper 


US$ 0.098 


1 91 9 Ink Cost per Newspaper 


US$ 0.010 


1 920 Total Cost (per Newspaper) 


US$ 0.108 


1 921 Cost per 1 000 Papers 


US$ 108.14 


1 922 Make Ready (per Week) 


23 


1 923 Savings per Week (Number of Papers) 


23000 


1 924 Savings per Week (Tons of Newsprint) 


4.5 


1 925 Savings per Year (Tons of Newsprint) 


235 


1 926 Savings per Week (US$ paper & ink only) 


US$ 2487 


1 927 Savings per Year (total US$) 


US$ 129,336 
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FIG. 30 





Image Based Pre-Registration Process 



FIG. 31 

Start 

Mark Recoginition 



> 



/ 



RGB Image Input 



7 



X 



Create the Binary Image (Bl) by 
thresholding with the 
"Starting Threshold" (ST) 
Bl= (R < ST) OR (G < ST) OR (B < ST) 



3100 
3101 

3102 



Error 



Create Linked List of possible marks: 
GOODLIST Creation 



T 



Find DUPLICATE Mark(s): 
Use Information from the GOODLIST 



Error 



Find Non-DUPLICATE Marks: 

Progressively raise threshold and 
scan for more marks inside the pattern 
tvindow circumscribing the Duplicate marks 



Error 



T 



Match Marks to Stations: 
Make a one-to-one correspondance 
between the detected marks and the 
printing stations defined in mark pattern. 



Error 
3109 

1_ 



T 



Check between Marks: 

Scan for any non-mark printed 
material between the marks. 



Exit 

Unsuccessful 



Exit 
Successful 
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FIG. 32 
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FIG. 41 
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FIG. 50 
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FIG. 51 



-3107 
FIG. 54 
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FIG. 32 



Start 

GOOD LIST Creation 
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3r 



I Calculate pixel dimension^ 



Get First Pixel P[n] from Bl at: 

x= WIDTH/3. y=0 



MARK COUNT = 0 



3206 



Yes 



Increase MARK 
COUNT by 1 




Process 
Vertical DIP 



3200 

3201 
■^F/G. 33 
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3204 
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FIG. 34 




FIG. 33 
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Start 

Calculate Pixel Dimensions 



3301 



Read XA' Pixel/Inch Constants (XP,YP): 

Camera Type selects lookup table. 
Zoom factor selects values 
in lookup table. 



3302 



3303 



3304 



3305 



Read WIDTH(ln.). LENGTH(in.): 

Obtain from Mark Specification. 



Calculate WIDTH, LENGTH: 

WIDTH =WIDTH(in.)*XP 
LENGTH= LENGTH(in.)*YP 



Read WIDTHVAR(ln.), LENGTHVAR(ln.): 

These are installation constants. 



Calculate WIDTHVAR, LENGTHVAR: 

WIDTHVAR= WIDTHVAR(in.)*XP 
LENGTHVAR= LENGTHVAR(in.)*YP 
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3307 



3308 



3309 



Calculate XMINDIP, YMINDiP: 

XMINDIP = WIDTH/10 
YMINDIP= LENGTH/10 



Read XBOUNDARY{ln.), YBOUNDARY{in.): 

These are installation constants. 



Calculate XBOUNDARY, YBOUNDARY: 

XBOUNDARY = XBOUNDARY(in.)*XP 
YBOUNDARY= YBOUNDARY(in.)*YP 
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Get next pixel on line below current: P[n] 



Exit 




3404 



Increment 
IVIARK 
COUNT 



3408 

DIP COUNT > YIVIIN DIP? 
3409 / ""^TTVes 

3411^ 

Does 

No MARK COUNT fall Yes 
inside LENGTH 
3410-—^ >>^Spec? 



3405 



Increment 

DIP 
COUNT 



Process 
Horizontal 
Direction 



MARK COUNT = 0 
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Exit 
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FIG. 35 

start 

Process Horizontal Direction 



Compute XA^ starting Coordinate: 

X = coordinate of vertical slice used in FIG. 32 
y = center of vertical section found in FIG. 32 







Find Riglit Edge 
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Find Left Edge 
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Compute width of mark: 
TESTWIDTH = RIGHT COUNT + LEFT COUNT 
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Find Right Edge 
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Increase RIGHT 
COUNT by 1 
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Process Right 
Horizontal DIP 
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Find Left Edge 
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Get First Pixel P[n] from Bl at: 

x.y 

4___^ 



LEFT COUNT = 0 
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Increase LEFT 
COUNT by 1 



Process Left 
Horizontal DIP 
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Start 

Perform Quick Boundary Check 




Calculate Bounding Box: 
From info from FIGs 30-35 
(xs.ys) (xe,ye) 



X 



Calculate Boundary Box: 

(x1 ,y1) = (xs-XBOUNDARY,ys-YBOUNDARY) 
(x2,y2) = (xe+XBOUNDARY,ye+YBOUNDARY) 



Get First Trixel Prgb[n] from RGB at: 

(x.y) = (x1,y1) 



X 



Save Trixel Prgb[n] in Srgb[n] 
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3802 

3803 
3804 



3805 



Get Next Trixel 
Prgb[n] from RGB: 

(x,y) = next pixel 
clock-wise around 
Boundary Box 



Add Bounding Box to GOOD BOX List 
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Process Right Horizontal DIP 




No 




RIGHT COUNT>0? 

Yes 



3903 



DIP COUNT = 1 



3900 

3901 
3902 



Yes 



More Right Pixels^ 
No 



3905-sr- Get next pixel to right of current: P[n] 



3904 



Increment 
RIGHT COUNT 
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Increment 
DIP COUNT 



Yes 





RIGHT COUNT = 0 
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Process Left Horizontal DIP 





LEFT COUNT > 0 ? 
Yes 



DIP COUNT = 1 



4000 



4001 



4002 



Yes 




4003 



4004 Get next pixel to left of current: P[n] 



4005 



Increment 
LEFT COUNT 




RIGHT COUNT = 0 
LEFT COUNT = 0 
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Exit 
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Start 

Find DUPLICATE Mark(s) 
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FIG. 42 
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Start 

Calculate Comparison Constants 



Read COLINEARITY(ln.): 

This Is an installation constant. 



4202 



Calculate XCOLINEARITY, YCOLINEARITY: 

XCOLINEARITY = COLINEARITY(in.)*XP 
YXCOLINEARITY=COLINEARITY(ln.)*YP 



4203^ 



Read COLORVAR[This is an installation constant]: 
This is the amount the DUPs interior color may vary. 



4205 



4204 



Read DUPX(ln.), DUPY(in.)[Obtain from Mark Specification]: 
These are the distances between duplicate marks, center-to-center. 



Calculate DUPX, DUPY, XSPECVAR, YSPECVAR: 

DUPX= DUPX(in.)*XP, DUPY= DUPY{in.)*YP 
XSPECVAR= (DUPX*SPECMATCH)/100, YSPECVAR= (DUPY*SPECMATCH)/100, 



Read WINDOWWIDTH, WINDOWLENGTH [Obtain from Mark Specification]: 
These are overall dimensions of a window which circumscribes the entire mark pattern 



4207- 



Calculate WINDOWWIDTH, WINDOWLENGTH: 

WINDOWWIDTH= WINDOWWIDTH(in.)*XP 
WINDOWLENGTH=WINDOWLENGTH(in.)*YP 



X 



4206 



Read XDUPBOUNDARY, YDUPBOUNDARY [These are installation constants]: 
This defines the area around the DUP mark(s) which must be clear of printed material 




Calculate XDUPBOUNDARY, YDUPBOUNDARY: 

XDUPBOUNDARY= XDUPBOUNDARY(in.)*XP 
YDUPBOUNDARY= YDUPBOUNDARY(in.)*YP 



4210, 



Exit 
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Start 
Single Column 



> 



Set Direction: 


s: 


VERTICAL 







Set DUP Boundary Checks: 

NORTH and SOUTH 



Find Colinear Marks: 

Given the set direction 
and checks, find two marks 
which work. 
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Start 
Single Row 



> 



Set Direction: 

HORIZONTAL 



X 



Set DUP Boundary Checks: 
EAST and WEST 



Find Colinear Marks: 

Given the set direction 
and checks, find two marks 
which work. 
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FIG. 48 
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Start 
Rows > Columns 



Set Direction and Test Box: 

VERTICAL, NO TEST BOX 

T 





Set DUP Boundary Checks: 




NORTHEAST and SOUTHEAST 


s 


OR 




NORTHWEST and SOUTHWEST 



X 



Find Colinear Marks: 

Given the set direction and checks, 
ind two marks which work (DUP1,DUP2) 



Set Direction and Test Box: 

HORIZONTAL. TEST B0X=DUP1 



Find Colinear Marks: 

Given the set direction and checks, 
find two marks which work (DUP3.DUP4; 



Set Direction and Test Box: 

HORIZONTAL, TEST B0X=DUP2 



I 



Find Colinear Marks: 

Given the set direction and checks, 
find two marks which work (DUP5.DUP6; 




J^^-'found DUP1^^6?> 



Exit 
Found 
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FIG. 48 



FIG. 48 
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Start 
Columns > Rows 



Set Direction and Test Box: 
HORIZONTAL, NO TEST BOX 



Set DUP Boundary Checks: 
NORTHWEST and NORTHEAST 
OR 

SOUTHWEST and SOUTHEAST 



Find Colinear Marks: 
Given the set direction and checks, 
Ind two marks which work (DUP1 ,DUP2) 



Set Direction and Test Box: 

VERTICAL, TEST B0X=DUP1 



Find Colinear Marks: 

Given the set direction and checks, 
find two marks which work (DUP3.DUP4; 



Set Direction and Test Box: 

VERTICAL, TEST B0X=DUP2 



Find Colinear Marks: 
Given the set direction and checks, 

Rnd two marks which work (DUP5.DUP6; 
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Get First Box on List: 

BOX 



Calculate BOX Center: 

(xcen.ycen) = BOX center 



Calculate TEST WINDOW: 

twx1 = xcen - WINDOWWIDTH/2 
twy1 = ycen - WINDOWLENGTH/2 

twx2 = xcen + WINDOWWIDTH/2 
twx2 = ycen + WINDOWLENGTH/2 



Get First Box on List: 
TESTBOX 



4705 




Yes 



Calculate TESTBOX Center: 

(txcen.tycen) = TESTBOX center 




4715 Yes 



Exit 
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4702 

4703 
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4708 



Get Next Box 
on List: 
TESTBOX 



DUP Boundary Checl( 

Check NORTH. EAST, 
SOUTH. & WEST of B0> 



Yes 
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Start 

Find Colinear Marks 




Get First Box on List: 

BOX 



4800 
4802 



Calculate BOX Parameters: 

(xcen.ycen) = BOX center 
(red, green, blue) = Average box color 
(width, length) = size of box 



Get First Box on List: 

TESTBOX 
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4806 4801 




4807 
4808 



Calculate TESTBOX Parameters: 

(txcen.tycen) = BOX center 
(tred,tgreen,tblue) = Average box color 
(twidth,tlength) = size of box 



Get Next Box on List: 

TESTBOX 



T 



Verify the following: 

ABS[red-tred] < COLORVAR 
ABS[green-tgreen] < COLORVAR 
ABS[blue-tblue] < COLORVAR 
ABS[width-twidth] < WIDTHVAR 
ABS[length-tlength] < LENGTHVAR 

If Direction=H0RIZONTAL: 

ABS[ycen-tycen] < YCOLINEARITY 
ABS[xcen-txcen] < XSPECVAR 

If Direction=VERTICAL: 

ABS[xcen-txcen] < XCOLINEARITY 
ABS[ycen-tycen] < YSPECVAR 



I 



4809 



Exit 
Found 



4814 




Good 
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DUP Boundary Check: 
Check set directions if no good try 
Alternate directions if set. 
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Start 

DUP Boundary Check 



4900 



Set Initial Boundaries: 

wBOUNDARY= WEST boundary 
eBOUNDARY= EAST boundary 
nBOUNDARY= NORTH boundary 
sBOUNDARY= SOUTH boundary 



4901 




Calculate Boundary Box: 
(x1.y1) = (xs-wBOUNDARY,ys-nBOUNDARY) 
(x2,y2) = (xe+eBOUNDARY,ye+sBOUNDARY) 



Yes 

Check (x2.y1) TO (x2,y2) | 





Check (x1.y1) TO (x2,y1) \ 



^est/SW/NW t 
Yes 



Check (x1.y1) TO (x1,y2) | 



905 




No 



Check (x1,y2)T0 (x2,y2) | 



4912 



Bump Out Boundaries 
that are enabled. 



4913 



Check Completed Boundaries: 

wBOUNDARY > XDUPBOUNDARY 
eBOUNDARY > XDUPBOUNDARY 
nBOUNDARY > YDUPBOUNDARY 
sBOUNDARY > YDUPBOUNDARY 




No 



4914 
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Start 

Find Non-DUPLICATE Marks 




Calculate the Region of Interest: 

(x1 ,y1 ) - (x2,y2) region which 
circumscribes DUP marks. 
WINDOWWIDTH, WINDOWLENGTH 



Set Starting Threshold: 

T = ST 



i 



Create the Binary image (Bl): 

Bl= (R < T) OR (G < T) OR (B < T) 



5004 



T 



Create Linked List of possible marks; 
GOODLIST Creation 



500 




Yes 



5006~-5:— 



No 



Bump up Threshold: 

T = T + BUIV1P 
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Exit 

Not Enough Marks 
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5001 



5002 



5003 



FIG. 32 
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Exit 

Enough Marks 
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Start 

Match Marks to Stations 



5100 



5104 



Calculate the Region of Interest: 

(x1 ,y1 ) - (x2,y2) region which 
circumscribes DUP marl<s. 
WINDOWWIDTH. WINDOWLENGTH 



Match Single Column 
or Rows > Columns 



f 

FIG. 53 





Match Single Row or 
Columns > Rows 



5102 
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FIG. 52 
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Yes , 


Sample 

Inside 
Mark Image 
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Sample 

Outside 
Mark Image 
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Exit 
Bad Match 



> 
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Exit 
Good Match 
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Match Single Row, Columns 

i 



> Rows 



Divide Region into Rows 



5202 



Set Next Row as Region of Interest 



5205 



Match a Mark to Station: 
From Left to Right 



5200 
5201 



5204 




More \ Yes 

Rows?, 



FIG. 53 
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Start 

Match Single Column, Rows > Columns 




Divide Region into Columns 



Set Next Column as Region of Interest 
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Match a Mark to Station: 

From top to bottom 



FIG. 54 

Start 

Check between Marks 



> 



5400 



Calculate the Region of Interest: 

(x1 ,y1 ) - (x2,y2) region which 
circumscribes DUP marks. 
WINDOWWIDTH, WINDOWLENGTH 



5404 




5401 



5402 



Check Single Column 
or Rows > Columns 



I 



FIG. 56 



5403 



Check Single Row or 
Columns > Rows 



No 




T 
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Exit 
Bad Check 



Exit 
Good Check 
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Check Single Row, Columns 



5506 



Yes 




> Rows i 



5500 



Divide Region into Rows 



5508 
5511 



Set Region of Interest: 

(rx1,ry1)-(rx2.ry2) Row 



X 



Get First Station 
Matched Mark in Row 
(ml x1 ,m1 y1 )-(m1x2.m1 y2) 



5501 
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5503 



No 




5504 



Get Next Station 
Matched Mark in Row 
(m2x1 ,m2y1 )-(m2x2,m2y2) 



X 



Check Region for Print: 

(m1x2,cy1)-(m2x1,cy2)) 



5505 



5507 




Print ? 



Yes 



Exit 
Good Check 



No 



Save Mark Location: 
(nn1x1,m1y1)-(m1x2,m1y2) = 
(m2x1 ,m2y1 )-(m2x2,m2y2) 



Exit 
Bad Check 
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Check Single Column, Rows 
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> Columns i 



Divide Region into Columns 



5600 
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Set Region of Interest: 

(cx1 ,cy1 )-(cx2,cy2) Column 



T 



Get First Station 
Matched IVIark in Column 
(m1x1,m1y1)-(m1x2,m1y2) 



5602 
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Yes 



5607 



Get Next Station 
Matched Mark in Column 

(m2x1 ,m2y1 )-(m2x2,m2y2) 



Check Region for Print: 

(cx1,m1y2)-(cx2,m2y1) 



Save Mark Location: 
(m1x1,m1y1)-(m1x2,m1y2) = 

(m2x1 ,m2y1 )-(m2x2,m2y2) 
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Yes 



5610 



\^ Bad Check 1 
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Exit 
Good Check 
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Web Press 
Motor Control 
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Operator 
Interface 
Display 
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Web 
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Mark 
Recognition 



Light IVIark 
Recognition 
(optional) 



Adjust 
Web 
Press 



5805 



Operator 
Interface 
Display 
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Web 
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Mark 

Recognition 



Light Mark 
Recognition 
(optional) 



Adjust 
Web 
Press 



Move Camera 
5906 ^ or 

Select New Camera 



5905 



Operator 
Interface 
Display 
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6001 



Start 
Light Marks 
Detection Program 



6002 




RGB Image, 
Number marks = numMarks, 
min/max Mark Width = min/maxWidth, 
min/max IVIark Heiglit = min/max HeiglitJ 
Input 



6003 



6004 




6005 



Segment RGB/HSI 
Images to Obtain 
Light Marks Image (LMI) 



^ FIG. 67 
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Image Process LMI: 
Determine Max, Min Thresholds 
for enhanced LMI plane 
(Rl) 



FIG. 62 
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Label and Measure all 
Light Marks In Rl: 
result is list of CBB's 
nMarks long: Stored 
Mark Table = SMT 



T 
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Exit 

Light Marks Detection Program 
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start 

Light Mark Segmenter 
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RGB/HSI Image 

Inputs 



/ 



Get Maximum of R,G,B 
Write to Max Plane 
(MP) 



Gray Level Erode 
MP Plane 



T 



Gray Level Dilate 
MP Plane 



T 



Low Pass MP Plane 

T 



E 



Determine Lowest Threshold 
Value = Dthresh from 

Histogram Lowest Maximum 

* 



Dup Mark Map (DMM) = 
Threshold MP Plane at Dthresh 



AND Intensity (I) image with 
DMM = Light Marks Image (LMI) 



X 



Return 
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Start 

Image Process For Light Marks 



'Light Marks Image 
LMI Input 



Enhance LMI Plane 









Determine Substrate 


6206 s 


Intensity Intensity Threshold 




= Sthresh 




Y 




Determine Dark Threshold 


6207 s 


Value and divide 




by 2 = Dthresh 



6208- 



Return 
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Gray Erode Enhanced 
LMI Image = El 











Low Pass Enhanced 
LMI Plane (El) = Result 
Image (Rl) 











Histogram Rl Plane 
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FIG. 63 

Start 

Light Mark Detector 



) 



6301 




Rl, Sthresh, Dthresh, 
Number of Marks = nMarks, 
min/max Mark Width = min/maxWidth 
min/max Mark Height = min/maxHeightj 
Input 
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6304 

6305 



Set Thresh = Dthresh 



Threshold Rl at Thresh = LMM 









Label LMM = LCM 




i 




Sort LCM Largest 


s 


to Smallest 
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FIG. 64 
D 
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Set Current Bounding Box (CBB) = 
Get Bounding Box of Largest Label 
Get CBB Width = CBBWidth 
Get CBB Height = CBBHeight 
Get CBB Center = CBBCenter 






< 1 
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CBBWidth > 
minWidth && 
CBBWidth < maxWidth && 
CBBHeight > minHeight && 
^BBHeight < maxHeighJ 
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Stored Mark Table 
SMT 
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Set Thresh = Sthresh 



Threshold Rl at Thresh = LMM 
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Label LMM = LCM 



Sort LCM Largest 
to Smallest 
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Set Current Bounding Box (CBB) = 
Get Bounding Box of Largest Label 
Get CBB Width = CBBWidth 
Get CBB Height = CBBHeight 
Get CBB Center = CBBCenter 
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'CBBWidth < maxWidth &&^ 
^BBHeight > minHeight &&^ 
CBBHeight < 
maxHeight 



No, 
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Get First RGB 

Pixel: R[n],G[n],B[n] 



Calculate Nomnalized RGB: 
Sum = R[n] + G[n] + B[n] 
R = R[n]/Sum 
G - G[n]/Sum 
B = B[n]/Sum 
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Calculate Intensity 
l[n] = R + G + B 



Calculate Saturation 
S[n] = 1 - ( (3 * min (R.G.B) )/ l[n]) 



Calculate Hue 
Numerator = ((R - G) + (R - B))/2 
Denominator = sqrt((R-G)*{R-G) + (R-B)*(G-B)) 
H[n] = arccos(Numerator/Denominator) 
if B > G, then H[n] = 2pi - H[n] 
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Pmin = minimum(P[n-4], 
P[n-3].P[n-2].P[n-1],P[n],P[n+1]. 

P[n+2],P[n+3],P[n+4]) 
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Find Largest Pixel Value 
Pmax = maximum(P[n-4], 
P[n-3].P[n-2].P[n-1].P[n].P[n+1], 
P[n+2].P[n+3],P]n+4]) 
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kernel values: 
Ptot = P[n-4]*K[n-4] + P[n-3]*K[n-3] 
P[n-2]*K[n-2] + P[n-1]*K[n-1] + 

P[n]*K[n] + P[n+irK[n+1] + 
P[n+2]*K[n+2] + P[n+3]*K[n+3] + 
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HISTOGRAM IM: 

Min = minlmum_value(IM) 
Max = maxlmum_value(IM) 



Shift IM: 
Set Every Pixel Value (P[n]) 
in IM: P[n] = P[n] - Min 
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Stretch IM: 
Set Every Pixel (P[n]) in IM 
P[n] = P[n]*(255/(Max - Min)) 
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